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Predicting dislocation climb: Classical modeling versus atomistic simulations 
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that an infinite straight dislocation of Burgers vector b 
and character 9 climbs at a velocity 
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The classical modeling of dislocation climb based on a continuous description of vacancy diffusion 
is compared to recent atomistic simulations of dislocation climb in body-centered cubic iron under 
vacancy supersaturation [Phys. Rev. Lett. 105, 095501 (2010)]. A quantitative agreement is 
obtained, showing the ability of the classical approach to describe dislocation climb. The analytical 
model is then used to extrapolate dislocation climb velocities to lower dislocation densities, in the 
range corresponding to experiments. This allows testing of the validity of the pure climb creep 
model proposed by Kabir et al. [Phys. Rev. Lett. 105, 095501 (2010)]. 

PACS numbers: 61.72.Lk, 62.20.Hg 

Dislocations can move out of their glide planes through 
the emission or absorption of point defects. Such a mech- 
anism, known as dislocation climb, has been modeled for 
more than half a century now using a continuum descrip- 
tion of matter and diffusion theory-^. One can thus find 
in most of the textbooks on dislocations^— analytical ex- 
pressions which give the dislocation climbing velocity as 
a function of the applied stress, the temperature, the 
point defect supersaturation, etc. Although this clas- 
sical description has been shown to reasonably explain 
experimental observations^—, a quantitative validation 
by direct comparison to experiments is out of reach. Dis- 
location climb is indeed rarely the single mechanism pro- 
ducing plastic strain and one has usually to deal with a 
complex dislocation microstructure. 

On the other hand, atomistic simulations can be used 
to study the pure climb either of an isolated dislocation 
or of a well controlled dislocation microstructure. Such 
simulations thus offer a natural way for a quantitative 
validation of the classical approach. In this Brief Report, 
we compare predictions of the classical dislocation climb 
modeler— to the results of atomistic simulations published 
by Kabir et al. in Ref. H3. 

Kabir et al. 14 performed atomistic simulations to study 
the climb of a mixed dislocation in body-centered cubic 
iron. Their simulations are based on a kinetic Monte 
Carlo algorithm which reproduces the diffusion of vacan- 
cies and their annihilation on a dislocation jog. The in- 
fluence of the dislocation on the vacancy migration bar- 
riers is fully taken into account thanks to an empirical 
potential which allows them to search for the minimum 
energy path joining all neighboring vacancy stable po- 
sitions (nudged elastic band method). In these simula- 
tions, the dislocations can be considered as being satu- 
rated with jogs because of the small dislocation length, 
and the climb is driven by a high vacancy supersatura- 
tion. 

When pipe diffusion is fast enough and a high con- 
centration of jogs is present on the dislocation, one can 
assume that vacancies are at equilibrium all along the 
dislocation linear—. Diffusion theory then predicts^— : 15 : 16 



The geometric factor 7/ = 2ir for an isolated dislocation. 
Dy is the vacancy diffusion coefficient in a dislocation- 
free crystal, Cy q the concentration of vacancies in equi- 
librium with the dislocation at a distance r c ~ b from the 
line, and Cy° the average vacancy concentration imposed 
at a distance Roc- This should correspond to half the av- 
erage distance between dislocations, i.e. Roc — l/(2 v /p5") 
if pd is the dislocation density. For a high vacancy su- 
persaturation, like in the atomistic simulations of Kabir 
et al. 14 , Cy q < Cy° and then 
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Equation (JT|) predicts that the climb velocity, once nor- 
malized by DyCy°, should only depend on pp. This 
agrees with the atomistic simulations of Ref. [lj: Results 
for all temperatures and vacancy supersaturations are 
well reproduced by Eq. (jTJ) . (Figs. [1] and H]). The two pa- 
rameters rj and r c appearing in this equation were used as 
fitting parameters, and the best quantitative agreement 
with atomistic simulations was obtained for the values 
r] = 12.8 and r c = 4.36. 

The geometric factor obtained from this fit (rj ~ 47r) is 
close to its 2-7T theoretical value. The slight difference may 
come from the fact that a periodic array of dislocation 
dipoles has been modeled in the atomistic simulations of 
Kabir et al^, whereas an isolated dislocation is assumed 
in Eq. ([!}. Burke and Nb*il also showed that the elas- 
tic interaction between the vacancy and the dislocations, 
which is neglected in our modeling approach, leads to a 
value for r\ slightly higher than 2-7T, thus in agreement 
with what we obtained. 

The conventions are usually to take the capture radius 
r c equal to the norm of the Burgers vector. Our fit of 
Eq. ([TJ leads to r c = 4.3 b, which actually agrees with 
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FIG. 1. Variation with the dislocation density pu of the climb 
velocity v c i normalized by DvCy for different temperatures 
and vacancy supersaturations. Symbols corresponds to atom- 
istic simulations 14 and lines to Eq. ([T]l. 



the definition of the core region (r < 4 b) that has been 
deduced from previous atomistic calculations^ on the 
same model system. 

One therefore sees that the two parameters rj and r c are 
physical parameters whose precise values are close to the- 
oretical ones. It is also important to notice that both pa- 
rameters do not theoretically depend on the dislocation 
density, the vacancy supersaturation, the applied stress, 
nor the temperature. This was checked and all atomistic 
simulations of Ref. [HI could be reproduced with a single 
set of parameters. 

This comparison with atomistic simulations shows that 
the classical modeling of dislocation climb leads to quan- 
titative predictions. Such a model, based on a continuous 
description of vacancy diffusion, does not explicitly take 
into account all atomic details of the vacancy diffusion 
close to the dislocation. It nevertheless manages to per- 
fectly reproduce results of atomistic simulations. Finally, 
all atomistic details on the vacancy migration are not 
so relevant to model dislocation climb and one can use 
classical mesoscopic approaches 3 -"— i 15 i 16 leading to simple 
analytical expressions like Eq. (p}. 

One can now use this expression of dislocation climb 
mobility to compare predictions of different creep models 
with experimental data, in particular the creep model de- 
veloped by Kabir et al^. The dislocation densities in the 
atomistic simulations^ (pd > 10 16 m -2 ) are much higher 
than the ones usually observed in creep experiments in 
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(Pt> J; 10 m ), As a consequence, one cannot 
directly use the dislocation climb velocities measured in 
atomistic simulations, but one needs to extrapolate them 
to lower dislocation densities. Kabir et alM- used power 
laws to perform such an extrapolation and concluded to 



the agreement of their model with creep experiments. As 
such power laws do not rely on any physical ground, it 
is worth checking if the same nice agreement can be ob- 
tained when the classical modeling of dislocation climb 
leading to Eq. ([T]) is used for this extrapolation. 

In the creep model, one usually uses a power law v c \ oc 
Pd™ to reproduce the variations of the climbing velocity 
with the dislocation density. Equation ([T]) leads then to 
an exponent m which depends on the density pd 
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This exponent can now be used in the creep model pro- 
posed by Kabir et al 1 ^ . This model assumes that all the 
plastic strain in creep is produced by climbing disloca- 
tions. Orowan law gives then the creep rate e — pi)bv c \. 
The dislocation density is fixed by the Taylor relation 
(Pd = (cr/aGb) 2 where a ~ 0.4 is an empirical constant 
and G is the shear modulus). The vacancy supersatura- 
tion, which varies linearly both with the applied stress a 
and the dislocation velocity, is given by Eq. (1) in Ref. US. 
Equation ^ combined with these assumptions leads then 
to a steady-state creep rate whose stress and tempera- 
ture dependences can be reproduced by the power law 
e = A(T n cxp (—Q/kT) with a stress-dependent exponent 
n = 3 + 4m: 
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Equation (j3|) allows applying the creep model of Kabir 
et al. to a stress range much lower than in atomistic sim- 
ulations. Stress in creep experiments in iron usually does 
not exceed 100 MPa^—. Equation ((3]) then predicts an 
exponent n smaller than 3.5, thus far from the values 
higher than 6 found experimentally 2 -^. One therefore 
sees that the creep model proposed by Kabir et al. can- 
not explain experimental stress exponents measured in 
iron. It naturally leads at low stress to an exponent close 
to 3, like any other creep model based on pure climb 25 . 
The high value for n obtained by Kabir et al.— corre- 
sponds to the much higher dislocation density of their 
atomistic simulations and cannot be directly compared 
to these experimental values. 

One probably needs to consider both dislocation glide 
and climb to obtain a stress exponent close to the exper- 
imental one. A creep model based on dislocation pure 
climb is only valid when dislocations cannot glide be- 
cause of some constraints, crystallographic constraints 
for instance as in hep metal o 26 ' 27 or quasicrystal a 13 ' 28 . 
Creep in alpha iron is far from this ideal case as there 
is nothing preventing the dislocations from gliding. As 
pointed out by Weertman in his review paper—, "almost 
all of the creep strain is produced by glide motion of dis- 
locations" . Although dislocation climb is the rate limit- 
ing process, dislocation glide strongly affects creep and 
cannot be ignored. Weertma n 25 ' 29 ! 30 showed for instance 
that the consideration of glide in a creep model makes 
the stress exponent increase from n — 3 to n = 4.5. 
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FIG. 2. Variation with the dislocation density po of the climb velocity v c i for different temperatures and for a vacancy 
concentration (a) Cy 3 = 510 -5 , (b) = 10~ 4 , and (c) Cv° = 10~ 3 . Symbols corresponds to atomistic simulations^ and 
lines to Eq. {!]). 



Another discrepancy between the creep model of Kabir 
et alJ^ and experiments^— in iron comes from the tem- 
perature dependence of the stress exponent n and the 
stress dependence of the creep activation energy Q. Ac- 
cording to Eq. ([3]), n does not depend directly on the 
temperature. A temperature dependence may only arise 
through a variation of the capture radius r c with the 
temperature^. As shown above, a single value of r c could 
be used to reproduce all the dislocation climbing veloc- 
ities obtained by atomistic simulations in the tempera- 
ture range 800 - HOOK. As a consequence, this creep 
model does not lead to any temperature dependence of 
the stress exponent n. The creep rate is obtained from 
Orowan law i = pubv c \ where the dislocation density is 
deduced from Taylor relation and the climbing velocity 
is given by Eq. ([lj . The activation energy Q for creep is 
then the activation energy for vacancy diffusion whatever 
the applied stress. The dependences found by Kabir et 
al. (Fig. 4 in Ref. ll4T ) are artifacts caused by a fit of the 
climbing velocity at high densities using a simple power 
law which neglects the logarithm appearing in Eq. ([1}. 
To explain the dependence observed experimentally^—, 



one has therefore to rely on mechanisms which are not 
considered in the creep model of Kabir et al.. 

Finally, it is worth pointing that the stress enters in 
the creep model of Kabir et al^ only through the control 
of the dislocation density and of the vacancy supersatu- 
ration. As the climbing velocity does not depend then 
on the dislocation orientation, all dislocations, whatever 
their orientations, are climbing at the same velocity. If 
there is no specific climbing direction being enhanced by 
the stress, no average macroscopic strain can develop and 
there will be no creep. 

In summary, the comparison with atomistic simu- 
lations shows that a classical approach at a meso- 
scopic scale manages to quantitatively describe disloca- 
tion climb. Such an approach not only allows rationaliz- 
ing results of atomistic simulations, but it is also neces- 
sary to extrapolate them in a range of dislocation den- 
sities corresponding to experiments. Thanks to such an 
extrapolation based on a physical sound model, a fair test 
of the validity of the creep model proposed by Kabir et 
alJ^ can be made, thus showing its inability to reproduce 
experimental data on creep in iron. 



* emmanuel.cl ouet@cea.fr 



J. Appl. Phys., 26, 1213 (1955) 



N. F. Mott, [Proc. Phys. Soc. B. 64. 729 (1951) 
J. Weertman 
J. Friedel, Dislocations 
J. P. Hirth and J 



(Per gamon Press, Oxford, 1964). 
Lothe, Theory of Dislocations, 2nd ed. 

( Wiley, New York, 1982). 

D Caillard and J. L. Martin, Thermally Activated Mecha- 



nisms in Crystal Plasticity, edited by R. W. Cahn (Perga- 
mon, Amsterdam, 2003). 

F. R. N. Nabarro, Theory of Crystal Dislocations (Oxford 



Fields and Dislocation Mobility, Modern Problems in Con- 
densed Matter Sciences, Vol. 31, edited by V. L. Indenbom 
and J. Lothe (North-Holland, Amsterdam, 1992) Chap. 10, 
pp. 571-624. 
J. 



Silcox and M. J. Whelan 
D. N. Seidman and R. W. BaUuffi 



Philo s. Mag.. 5. 1 (1960) 



531 (1966) 



Phys. Status Solidi, 17. 



uffi, [Phys. Status Solidi B, 31, 443 (19691 



R. W. Bal 

J. Powell and J. Burke, [Philos. Mag., 31, 943 (1975)| 
C. Bonafos, D. Mathiot, and A. Claverie, J. Appl. Phys. 



U liv. Press, London, 1967). 
V. L. Indenbom and Z. K. Saralidze, in Elastic Strain 



83, 3008 (1998) 



F. Mompiou and D. Caillard, Acta Mater., 56, 2262 (2008) 



4 



M. Kabir T. T. Lau. D. Rodney. S. Yip, a nd K. J 



Van Vliet, [Phys. Rev. Lett., 105, 095501 (2010) 
15 D Mordehai, E. Clou et, M. Fivel, and M. Verdier, iPhilos 



P. W. Davies. G. Nelmes. K. R. Williams, and B. Wilshire, 



M |ag., 88, 899 (2008) 
16 B Bako. E. Clouet. L M. Dupuy, and M. Bletry, iPhilos 



Met. Sci. J., 7, 87 (1973) 



M |ag., 91, 3173 (2011)1 
17 M. A Burke and W. D. Nix, |Philos. Mag. A, 37, 479| 



(1978) 



J. Weertman, in Rate Processes in Plastic Deformation of 
Metals, edited by J. C. M. Li and A. K. Mukherjee (Am. 
Soc. Metals, Metals Park, Ohio. 1975) pp. 315-336. 
G. Edelin and J. P. Poirier, |Philos. Mag., 28, 1203 (1973' 



R. Le Hazif, G. Edelin, a nd J. M. Dupouy, 



Mater., 60, 399 (2009) 



Trans. B, 4, 1275 (1973) 



Metall. Mater 



T T. Lau. X. Lin. S. Yip, and K. J. Van Vliet. |Scripta| za F. Mompio u and D. Caillard. pVlater. Sci. Eng. A, 483-484 



F. Garofalo, L. Zwell, A. Keh, and S. Weissmann, Actaf 



~7T 



143 (2008) 



and H. Oikawa, 



M |etall., 9, 721 (1961) . 
S. Karashima, T. Iikubo, T. Watanabe, 
Trans. Jpn. Inst. Met., 12, 369 (1971). 
S. Karashima, H. Oikawa, and T. Watanabe, Acta Metall 
l 4, 791 (1966)| 



J. Weertman, ASM Transactions Quarterly, 
(1968), ISSN 0097-3912. 
J. Weertman, 



61, 681 



Phys. Rev., 107, 1259 (1957) 



J. Cadek and K. Milicka, |Czech. J. Phys., 18, 1156 (1968' 



J. Cadek. M. Pahutova . K. Ci'ha, and T. Hostinsky, 
Mjetall., 17, 803 (1969). 



Acta 



The temperature dependence of the shear modulus G is 
taken into account by normalizing the applied stress a with 
G when processing results of creep experiments. A tem- 
perature dependence of the parameter a entering Taylor 
relation is also possible but the modeling of such a varia- 
tion will require a more complex approach than the simple 
creep model used here. 



